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TECHNICAL NOTE 3283 


AERODYNAMIC FORCES, MOMENTS, AND STABILITY 
DERIVATIVES FOR SLENDER BODIES OF 
GENERAL CROSS SECTION 
By Alvin H. Sacks 

SUMMARY 


TLe problem of determining the total forces, moments, and stability 
derivatives for a slender body performing slow' maneuvers in a compressible 
fluid is treated within the assumptions of slender -body theory. General 
expressions for the total forces (except drag) and moments are developed 
in terms of the geometry and motions of the airplane, and formulas for 
the stability derivatives are derived in terms of the mapping functions 
of the cross sections. 

All components of the motion are treated simultaneously and second 
derivatives as well as first are obtained, with respect to both the 
motion ccmponents and their time rates of change. Coupling of the longi- 
tudinal and lateral motions is th\is automatically included. A number of 
general relationships among the various stability derivatives are fovind 
which are independent of the configuration, so that, at most, only 35 
of a total of 325 first and second derivatives need be calculated 
directly. Calculations of stability derivatives are carried out for two 
triangular wings with camber and thickness, one with a blunt trailing 
edge, and for two wing-body combinations, one having a plane wing and 
vertical -fin. 

The influence on the stability derivatives of the squared terms in 
the pressure relation is demonstrated, and the apparent mass concept as 
applied to slender-body theory is discussed at some length in the light 
of the present analysis. It is shown that the stability derivatives can 
be calculated by apparent mass although the general expressions for the 
total forces and moments Involve additional terms. 


INTRODUCTION 


Ever since R. T. Jones (ref. l) in 19 ^ demonstrated the use of 
Munk^s apparent mass concept of 1924 (ref. 2 ) for solving problems of 
slender wings in a compressible flow, an ever-increasing number of 
investigators have entered the field of analysis now commonly known as 
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slender-body theory. The stability derivatives of slender triangular 
wings were treated by Ribner (ref. 3) in 19^7 following the pattern of 
Jones, and in 1948 Sprelter (ref. 4) extended the latter's result by 
means of conformal mapping to include certain wing-body combinations. 
Shortly thereafter, in 1949, Ward's general analysis for slender pointed 
bodies in steady supersonic flow (ref. 5) was published. 

After. the appearance of Ward's analysis, a number of papers were 
written on various aspects of slender-body theory including extensions 
to subsonic flow and to "not -so -slender" bodies (e.g., refs. 6 and 7), 
and in 1952 Phythian (ref. 8) develojped an analysis which included time 
variations in forward velocity and angles of incidence. Although many 
papers (e.g., refs. 9 and lo) have been devoted to the calculation of 
various stability derivatives for specific configurations, it is only in 
the past few months that a report by Miles (ref. 11) has given the com- 
plete counterpart of Ward's analysis for unsteady flow. 

The determination of stability derivatives has long been of concern 
to the engineer in connection with the dynamic behavior of airplanes, 
but the problem has assumed even greater proportions in the more recent 
slender configurations of missile design. The stability derivatives them- 
selves correspond to the coefficients of a Taylor expansion representing 
a particular component of force ( say lift) or. moment as a function of 
the airplane motions. The coefficient of any particular motion (say q) 
in the expansion is equal to the partial derivative of the force or 
moment component with respect to that motion. Ordinarily, stability 
derivatives are defined as these partial derivatives evaluated with all 
of the independent variables except a set to zero, so that the usual 
stability derivatives depend upon the initial angle of attack as well as 
on the configuration. In the present paper, however, all derivatives 
^e^ evaluated with all of the independent variables (a, 3, p, g, r, d, 3, 
p, i") set to zero. The advantages of this choice will become apparent 
in the course of the analysis. 

The present paper employs an approach believed to be novel in 
slender-body theory and is concerned with developing formulas for the 
forces and moments as well as the stability derivatives for general 
slender wing -body combinations.^ The significance of the squared terms 
in the pressure relation for slender configurations precludes the 
possibility of considering the longitudinal and lateral motions inde- 
pendently, so all motions of the airplane are. treated simultaneously. 


■H/hile the present analysis was being carried out, Bryson (ref. 12) 
published a paper treating essentially the same subject from a different 
viewpoint based on the tacit assumption that all the forces, moments, 
and stability derivatives can be obtained from the apparent mass analogue. 
This assumption and some of Bryson's results are discussed in a later 
section. 
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The mathematical restrictions on the generality of the shapes that 
can be rigorously handled have been discussed in detail by Ward (ref. 5} 
and more recently for the unsteady case by Miles (ref. ll) . Such dis- 
cussion will not be repeated in this report. 
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LIST OF IMPORTAIJT SYMBOLS 


coefficient of — tenn in series expansion of the mapping 

00 


function £ = f(n) “ ^ 


n=o 


coefficient of p term in expansion of the complex potential 

f(£) ^ 

value of Ai at a= 3 =p = q^ = r = 0 
coefficient of Zn £ in expansion of 


distance from airplane nose to pivot -point 
complex potential cp + i'^ 
length of airplane 

force in the z direction (approximately lift) 
rolling moment about the x axis 
reference length 

pitching moment about pivot point x = ci 
yavlng moment about pivot point x = ci 
angular rolling velocity about the x axis 
pressure 

angular pitching velocity about the y 'axis 

fluid speed relative to axes fixed in the body 

component of normal to body contour in plane x = const, 

(positive into the fluid) 
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if 

component of tangential to body contour in plane 

X = const, (positive counterclockwise looking upstream) 

R V + IW 

r angular yawing velocity about the z" axis 

ro radius of transformed circle corresponding to airplane cross 

section 

S cross-sectional area 

reference area 
t time 

Uq component of flight velocity along negative x axis 

Vq component of flight velocity along positive y axis 

V Vq - r(x-Ci) 

Vi speed of a point, fixed in the xyz system of axes 

Ur^Vr^Wj. components of in the x,y, z directions 

Wq component of flight velocity along positive z axis 

W Wo - q(x-ci) 

V force in the y direction (side force) 

xyz Cartesian coordinates fixed in the body (x rearward, 

y to starboard, z upward) 

a angle of attack (angle between arbitrarily chosen xy plane 

and flight direction) 

3 angle of sideslip (angle between xz plane and flight 

direction) 

e angle between the positive y axis and the tangent to the 

body contoTor in a plane x = const. 

p fluid mass density 

V outward normal to the body contour in plane x = const. 

£ 


y + iz 
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£(, complex coordinate of centroid of cross-sectional area 

(yc + iZ(j) 

a complex coordinate in transformed circle plane 

9 velocity potential 


'Px> V 
^z^'Pt 


partial derivatives of 9 with respect to x,y,z, and t 




velocity potentials for unit velocity of the cross section in 
the z directions 


q>’ velocity potential associated ^Ith variations in shape and size 

of cross section with x 

i|f stream function 

stream function along the contour of the cross section 

5 


Special Notations 


/ 


contour integral taken once round the body cross section in the 
positive (cauntercloc3swise) sense 


Force coefficients; Cy = 


Moment coefficients; Cm = 


(l/2)pUo^Sr 
M 


(i/2)pUq^Si»Zj. 


> etc . 

etc • 


Stability derivatives: 
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R real part 

I imaginary part 

( ) derivative of ( ) with respect to time 

( ) complex conjugate of ( ) 

GENERAL AI^ALYSIS 


The problem to be treated here is the determination of the aero- 
dynamic forces and moments (except drag) and the stability derivatives 
for a smooth slender airplane or missile of arbitrary cross section per- 
forming slow maneuvers in a compressible fluid. The configuration will 
be limited in that the base (if any) of the fuselage and all wing trail- 
ing edges must lie in a plane essentially normal to the longitudinal body 
axis. : - 


Differential Eq.uation and Pressure Relation 


The linearized differential equation for the velocity potential of 
unsteady motion of a ccanpressible fluid is the well-known wave equation 

®TT ■ ° 

where the system of axes Ajig is fixed relative to the undisturbed fluid, 
Cq is the speed of sound in the undisturbed fluid, and t is time as 
measured in the system. Thus the velocity potential $ is express- 

ible as ^ = *^A,)i,g,T). 

In general, the pressure relation associated with the velocity 
potential $ is given by (ref. I3, p. 19) 

^ - ®T - I + f (2) 

where is pressure and q^^ is the magnitude of the fluid velocity 

expressible as^ 

• ( 3 ) 


^The subscripts on p and q are used to, distinguish them from the 
angular velocities of rolling and pitching. 
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It will be convenient for 
ttie present problem to introduce 
a coordinate system x^yz which 
is fixed in the airplane. The 
axes chosen for this pui*pose- are 
shown in the sketch and comprise 
a Cartesian coordinate system 
endowed with the translational 
velocities Uq, Vo. Wq and the 
rotational velocities p; q.^ r 
of the airplane. (Note that 
this does not constitute a com- 
pletely ri^t-hand system.) 

The Xi axis passes through 
the airplane nose, and the ori- 
gin of the xiyz system is 
fixed at an arbitrary distance 
Cl from the nose as shown in 
the sketch. 





Since it is the purpose of this paper to study only instantaneous 
forces and moments (l.e., no time histories), it will be sufficient to 
choose an instant of time such that the positions of the moving Xiyz 
system and the stationary A|J.^ system are just coincident. Thus, 
equations (l) and ( 2 ) will be expressed in the xiyz system only for 
this instant, designated r = 0. For this purpose a new function T 
is Introduced such that 


'PCxijyjZjt) = 0(>v,|a,|,T) (4) 

Now, through the use of the transformations relating the moving and 
stationary coordinates (see e.g., ref. I 3 , p. 12 , and ref. l4, p. 79) 
one finds at t = t = 0 that 

^ ^ ^ ^ ^xi ^ ^9 ^y j 

^ ^t ^Xi Sy 'dr 

= + UqTxi *" + pz - rxi)cpy - (Wq - py - qxi)9z 


and 


= = 'Py; = 

It can be seen from the sketch that the quantities (Vq + pz - rxi) and 
(Wq - py - q.xi) are simply the velocity components in the y and z 
directions of a point fixed in the xiyz system. Note that in the 
corresponding x component (- Uq + ry + qz) the products ry and qz 
are considered negligible compared with Uq* 
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With the assumptions of slender bodies, small singles, and slow 
maneuvers, the differential equation (l) reduces to Laplace* s equation 
in planes xi = const, near the body; that is 

^yy ^ZZ ~ ^ 

It follows also that the density p must be treated as a constant 
in the pressure relation (2) which now becomes 






+ $1^) + const. 


= - ^t - ^o*^Xi + ('^o + pz - rxi) 9 y + (Wq - py - qxi) 9 z “ 

I ('Py® + 9 /) + const. (6) 


Here, for the sake of convention, one can transfer the origin of the 
moving axes to the body nose by letting 

Xi = X - Cl 

so that pitching and yawing rotations are still made about an arbitrary 
pivot point X = Cl. Thus, introducing the notation 

V = Vq - rxi = Vq - r(x - Ci) 

W = Wo - 4x1 = Wo - q,(x - cx) 
equation (6) can be written 

^ = - 9^ - Uq 9 x + (V + pz) 9 y + (W - py)?^ - i ( 9 y^ + + const. 

( 7 ) 

This, then, is the pressure relation (referred to the moving body axes) 
upon which the calculations of the forces and moments will be based. It 
should be noted that a consistent application of the slenderness approx- 
imation requires the retention of the squared terms and 9 z^. Thus 

slender-body theory is not a strictly linear theory although the differ- 
ential equation (5) is certainly linear. This means that solutions of 
equation (5) for cp (and hence the velocities) can be obtained by super- 
position, but the pressures cannot. Likewise, the forces and moments 
cannot be calculated by superposition except for those special cases in 
which the contribution of the squared tenns to the loading vanishes. 
Furthermore, when the airplane is performing combined maneuvers (e.g., 
simultaneous rolling and pitching), the squared terms may contribute 
additional forces and moments. These in fact give rise to the second- 
order stability derivatives that will be included in the present analysis. 
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Total Forces and Moments 


The analysis to be presented here 
for calculating the total forces 
(except drag) and moments on a slender 
configuration will be the counterpart 
of a method originally due to H, Blasius 
(ref. 15 ) for obtaining the forces and 
mcMents on a two-dimensional body of 
arbitrary shape immersed in a steady 
incompressible stream.® This analysis, 
although not suited to the calculation 
of total drag, will nevertheless take 
proper account of the local forces 
associated with leading-edge suction. 
Consider a lamina of the slender air- 
plane, cut parallel to the yz plane, 
of thickness dx as shown in the sketch. 
One can write immediately the differ- 
ential lift and side force on an ele- 
mental area in terms of the local 
pressure pi on the body: 


z 



and 


d^L = pi dy dx 


( 8 ) 


d^Y = - pi dz dx j 

Now, by introducing the ccxnplex variable ^ = y + iz, one can express 
the differential complex force as 

d^Y - i d^L = - Pi dz dx - i Pi dy dx = - i pi dx df (9) 


where ^ is the complex conjugate of ^ . In a similar fashion, the 
differential rolling moment about the x axis can be expressed as 

d^L’ = - Pi z dz dx - Pi y dy dx = - Pi dx R(^d5) (lO) 


where R denotes the real part. Further, the differential yawing and 
pitching moments about the pivot point x = ci are given by 

d?N - i d^M = - (d^Y - i d^L) (x - Ci) = i Pi(x - Ci) dx dl (ll) 

Integration of equations (9), (lO), and (ll) gives for the total forces 
and moments 


®The method of Blasius has been extended to two-dimensional unsteady 
incompressible flows by L. M. Milne-Thomson (ref. I6) and recourse will 
be had to many of his techniques throu^out the present analysis. 
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where the contour integrals are taken round the boundary of the airplane 
cross section In the positive (counterclockwise) sense. Before these 
integrations can be effected, the pressure pi must, of course, be 
expressed as a function of the complex variable Toward this end, 

it will be convenient to Introduce two new definitions pertaining to 
velocities in the fixed and moving coordinate systems. First, the square 
of the speed of a point fixed in the xyz system can be written as 

Vi^ = + (V + pz)^ + (W - py)® (13) 

Second, it is noted that the square of the fluid speed relative to the 
xyz system is given by 

q.r^ = Ur® + Vj? + Wr® = ( + Uq)® + (tPy - V - pz)^ + (tp^ - W + py)^ 

(l*t) 

SO that, neglecting 9 ^^ in comparison with <Py^ and 9 ^^^, one can write 
equation (?) in the form 


^ \ - I 1r^ + I ■'^1^ + const. (15) 

This expression will now be formed as a function of £ through the 
introduction of the complex variable R = V + iW and the complex 
potential F = 9 + 1 +. 


The speed Vi 
= Uo^ 

= Uo^ 


is immediately expressible as 


(v + pz) + i(w - py)l [”(v + pz) - i(W - 


J L 


py) 


+ (r - ip£) (R + ipD 


( 16 ) 


while the components v^ and Wj. of the relative fluid velocity are 
related to F through the complex velocity by 


vr 


dF — •» 

I’M'r = - (V + pz) + 1 (W - py) = ^ - R - 1 p 5 


(17) 



NACA IN 3283 


11 


Furthermore, at the body sur- 
face it will be seen from the 
sketch that 

+ iwp = + 


<lne 

= (<ls - 

so that 

Vr - iwj, = (qg + 

( 18 ) 

where and are the 
tangential and normal com- 
ponents of the transverse 

relative velocity and 0 Is the angle defined. in the sketch. Com- 
parison of equations (IT) and (18) gives 



(19) 


from which 

( S " ® ^ 

It is now noted that (see sketch above) 

■<■ 4/ = - Uo^ - 2Uo<Px 
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Finally, the boundary condition 
at the body siirface requires that the 
normal component of the fluid velocity 
relative to the moving xyz system be 
equal to the normal component of the 
forward flight velocity. Let C be 
the contour of the airplane cross 
section in .the plane x = const, and 
let Cl represent the projection on 
that plane of the contoiir at x + dx 
( see sketch) . If V is the outward 
normal to C at any point, and dv 
is the distance between C and Ci 
measured along the normal, then the 
above boundary condition is, to the 
present order of accuracy, 

= Uo g (21) 

Thus the pressure relation of equation (15) can now be written (for 
points on the airplane surface) in the desired form 

i (R - lp 5 ) (R + ipt) + const. ( 22 ) 

Reduction of the Integrals . - Before making' use of equation ( 22 ) in 
writing the integrals for the forces and moments, it will be useful to 
notice from the sketch on page 11 that the differential distances on the 
body contour in planes x = const, are related by the angle 9 so that 
dy = ds cos 9 and dz = ds sin 6 where ds is the differential arc 
length, positive counterclockwise. Hence, 

d£ = ds ei0; dl = ds e-10j dl = d£ e-2l0 (23) 

so that the first integral of equation (l2) for the complex lateral force 
can be written, after expanding the squared term in the pressure relation, 
in the form 


z 
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Y - IL = ip <P^d£ + IpUo dx f cp^j^df + 


J 

1 1 yj'dx f (g)"d£ i i I dx / (i i ipif d£ . 

ip J' RdxjT^d^+pp J' dxjTr^a£ + 
pUo y^^ ^ ^ ^ J'^ ^ f ^ * ip£)ds - 

i ^ jT dx jT (R - IpS) (r + ip 5 )dT { 2 k) 


Note that the constant term in the pressiire relation contributes nothing 
to the contour integrals of equation (12). 

The nine contour integrals of equation ( 24 ) can be divided into 
three types: (a) integrals that do not depend on the velocity potential; 

(b) integrals containing the real potential T; and (c) integrals con- 
taining the complex potential F. The first type can be integrated at 
once and these -will be dealt with first. The second type will be reduced 
to Integrals of the third type by determining the stream function on the 
boundary, and the third type will then be handled by the method of 
residues. 


It is first noted that P — ds is simply the rate of change of 

J dx 

/ — dv 

t — ds is the complex conjugate 
dx 

of the rate of change (in the x direction) of the moment of cross- 
sectional area. Thus, one can write 

*lpi(BT,) ( 25 ) 

where ^ is the complex conjugate of the position of the centroid of 
area of the cross section. 

The other two integrals of equation ( 24 ) that do not depend on the 
velocity potential can be conveniently evaluated by the use of Stokes’ 
theorem which can be stated in complex form as (see ref. 16, p. I30) 
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or 



(26) 


■where S is the area enclosed by the contour. Thus;, using the first 
form of equation (26), one can -write 


jT (R + IpD^d? = 21 



2ip(R + lpl)dS = - 4pS(R + ip^^) 


(27) 


where S is the airplane cross-sectional area. Similarly, from the 
second form of equation (26) , 


jT (R - ipt) (R + lp5)d£ = - 2pS(R + 1 p5(.) (28) 

Thus, all the Integrals of the first type discussed above have been 
evaluated. 

Before introducing the stream function for the evaluation of the 
first two integrals of equation (2i^-), it is well to note that the time 
differentiation can be taken outside the Integral sign with no diffi- 
culty, but the X differentiation cannot since the contour of inte- 
gration is Itself a function of x. It is shown in Appendix A that 


where C is the contour of integration round the airplane cross section 
eind the surface of the airplane can be expressed either as 

z = Zi(x,y) 


or 


y = 7i(x,z) 

Corresponding to these expressions for the surface are the expressions 
for the slopes of the surface 


^ dv/dx 
^x cos 6 

^Yi _ dv/dx 
sin e 


and 
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SO that, recalling the relations 



of equation (2l) . That is, the 

total outward nonnal of the fluid velocity in the plane x = const, is 
given by the sum of and the normal velocity of a point on the 

boundary considered fixed in the xyz system. Hence, (see sketch), 

= (V + pz)sin e - (W - py)cos e + Uq — (30) 

OS dx 

and it is recalled that sin G = dz/ds and cos 0 = dy/ds. The sense 
of ds is indicated by the arraw along the contour. Thus, integrating 
along the contour, one finds that the stream function on the surface is 
given by 

tg = Vz + I* pz^ - Wy + i 4- Uo J" ^ ds + G(x,t) 

where G(x,t) is an arbitrary function of x emd t. Now can also 

be expressed as a function of the complex variables £ and R by noting 
that 


KR^) = Vz. - Wy = - I (r5 - R^) 


and 


+ 2^ = U 
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whence 


% = - k rC) + r + G(x,t) 


1 _ 


dv 


dx 


(31) 


The integral ^ TdJ is now expressible in t_erms of the complex poten- 
tial since T = F + it and the integral jT td^ can be evaluated from 

equation ( 3 l)« That is, 'if one sets t = to + ti where ti represents 
the integral term of equation (31)^ then 


/ ■ 


ta£ = J' ilfocit; + jT tidS = jT - jT (32) 

Now^ ti taken once round the contour has the value Uq (f — 6 lb or 

J (ix 


dx 


dx 


dx 


simply Uq ^ and it is recalled that <j> C that 

upon evaluating ^ tod^. by the second form^of Stokes' theorem (e(i, (26)) 
and noting that jT G(x,t)dt = 0 , one finds 


/ 


>|fd£ = -S(R + ip£(.) + Uo^o — - Uo — (sEc) 

dx dx 


( 33 ) 


where Sq conplex conjugate of the point at which the integra- 

tion was begun on the boundary. The final expression for jT cpd^ 
is therefore 


9 d 5 =jT 5 d 5 '+ 1 jT yd^ =jTFde+ iUo^o^ - iS(R + iplg) - iUo-|^(s£o) 

(3l») 


and the time derivative is 
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Of the remaining integrals in equation (24) there is one which still 
requires modification before the method of residues can be applied. 
Specifically, since F = F + 2it one can write 


J' iSF = J' + 2idi|r) = jT + 21 jT 


Sat 


and the second integral can be integrated by parts using equation ( 33 ) 

and noting that ijr taken once round the contour has the value Uq — - 

dx 

The resulting expression is 

jT Sap = jT Sdi + 21 S(K + ip?^) + 2iUo ^ (sT^) ( 36 ) 


so that equation (24) for the complex force can finally be written 
(after collection and cancellation of terms) in the form 


Y - 


IL = ipy^ dx FdS + PjT S ^ dx + IpUo jT Fdl - 

ipUo jT id? - pUc 

PP J' dxj^ SdF,+ ipp J' S(R + ipSQ)dx ( 37 ) 


x=Z 

f - P^o TuoSo ^ - S(R + 2ipSg) - Uq (SSc) 

x=o L ^ dx 


It should be pointed out that several terms have vanished by virtue of 
the fact that the x axis passes throu^ the airplane nose. In par- 
ticular, note that S = = 0. at x = 0 . 

The complex potential .- Although equation (37) appears quite 
unwieldy, all of the contour integrals are now in a form which admits of 
evaluation by the method of residues. For a body moving through still 
air, as in the present problem, all velocities vanish at infinity and the 
complex potential F can be expanded in a Laurent series of the form 


F = B(x)ZnS + 


00 

Z 

n=i 




+ D(x,t) 


(38) 
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where B(x) is a source strength and the coefficients An(x,t) give the 
intensities of all the hi^er order singularities representing the 
desired body shape and motions. This expansion applies only for large 
values of but since there are no singularities outside the body, the 
contour integrals can be evaluated around a contour sufficiently large 
to Insure the validity of the expansion. The arbitrary function P(x,t) 
is of no concern here since it can contribute nothing to the contour 
Integrations. For the determination of drag, on the other hand, this 
function would be required. 

From equation .( 38 ) j the derivative of F is 


so that 



- S. 
d? " e 


^ s V 
~ i L 

n=i 



Now it is seen that the residue of (dP/dS)^ is zero since there is no 
l/t term in the expansion. Therefore, 


Also 

jT dF = y' ^ d£ = 2iti = 2i«lB 


(39) 


(»tO) 


Similarly, if one writes the conjugate function 

00 _ 


F=B(x)lnI+y + D(x,t) 

L. 5“ 

n=i 


it follows that 


and 


jT Fdl = - 2itiB5o - 2*iAi 


(^1) 


(42) 


+ 2itlAi 


(43) 
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The coefficient B (and, hence, B) can be evaluated by calculating the 

r 89 

integral 6 — ds in two different ways. Thus 
J on 

Mo E f 

and, by virtue of Gauss* theorem, 

21T 2ir 

/'|S^ds= f ride - f — ride = 2«B 

J dn Jq Jo ri 


where ri is the radius of a -circle enclosing the body cross section. 
Hence, it is seen that 


B =B =Ho3S 

23 t dx 


w 


The final exi)ression for the complex lateral force is obtained by 
using equations ( 39 ) to (44) in rewriting equation (37)- That is. 


r 


Y - IL = 2jtpUoA.i^^j + pUo S(R + SlpI^) + Uq ^ (S£o) + 


>£c>' 
-lx=l 


2jtp f dx. + p r S^dx + 2«lpp r Aidx + 

Jn 8t jU Ot UU 


a ^ 


_ 

ipp / S(R + ip?c)^ 




For the case of steady straight flight f ^ = p = 0 ) , the complex 

V 8 t 8 t 


force of equation (4^) reduces to that given by Ward (ref. 5)- Although 
equation (4^) applies to slender airplanes having cross sections of 
arbitrary shape, it is of interest ttiat in a large number of practical 
cases, it is possible to choose the x axis so as to place the center 
of cross-sectional area always along the axis and thus to make 
equal to zero. The simplest example would be an airplane having mirror 
symmetry of area about both the y and z axes. If the wings have no 
thickness, this places no restriction on the wings themselves with regard 
to number of wings, arrangement, dihedral, camber, etc. Equation (4^) 
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as It stands gives the total side force and lift on a slender configura- 
tion having a nose at the forward end. However, it can also be used to 
give the contribution to these forces of a segment of the airplane lying 
between the planes x = and x = by simply evaluating Y - iL for “ 

Z = Zi and Z = Zg and subtracting the results. 

If one considers the transformation of the arbitrary cross section 
in the £ plane to a circle in the a plane 

00 

I = f(cr) = a + ^ ^ (»t6) 

n=o 

it can be sho;^ (see Appendix B) that the dependence of the residue Ai 
(and, hence, Ai) on the rolling velocity p is determined by the form 
of f(a) and therefore by the shape of the airplane cross section. It 
is found that Ai is Independent of p if aQ vanishes and if f(a) 
contains either only odd negative powers of a (n odd) or else only even 
negative powers of d (n even). It can be seen that if n is odd then 
f(-a) = - f(cr) and, hence, the cross section has symmetry about two 
orthogonal axes. The statement can therefore be made that for airplanes 
having symmetry about both the y and z axes (no dihedral) one can 
determine the total complex force if he Imows only the complex potential 
due to pure translation in the yz' plane. 

An expression for the pitching and yawing moments can now be 
obtained -from equation (l2) by a procedure exactly parallel to that used 
in obtaining equation (4^) . Making use of th^ foregoing evaluations of 
the required contour ^ integrals , one finds that the resulting expression 
is . 


N - iM = - 2irpUo 


nl 

I 


Cl) 





S(R + 2iptc) 


+ 


Uo 


dx 




(x - Ci) ^ dx 
ot 



- Ci)s ^ ax - 


2itlpp r (x - ci)Aidx - Ipp r (x - Cl) (R + 1pSq)S dx (4?) 

Jq .. yo - . . 

The evaluation of the integrals for the rolling moment L’ is 
somewhat different due to the additional £ appearing in the expression 
of equation (12), so that the integrand of some of the contour integrals 
appearing here will be nonanalytic in the variable of integration. This 
precludes a direct application of the method of residues'. Such an 
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integral, which does, in fact, arise is I j! £ 5 dF. However, in this 
instance it can be shown with the aid of eguations (26) and (3I) that 

/ Kg a. (1.8) 

Evaluation of other required contour integrals by means of equations ( 26 ), 
(31), and (38) yields 





m 


jT SdF = - 2JtiAi 

and it can be shown in a manner parallel to that of Appendix. A that 

R/,.£ar-lt(|;/ 9 !d?)-l/g|U. ( 50 ) 

so that the final expression for the rolling moment can be written in 
the form 

L' = pUqR jT 9 ?d£ + pR dx ^ jT 9 td£ - 23TpI RAidx - 
x=Z 

ppR S(Rle - - PUoI ^ S (51) 


or, since R(£d0 =id(££). 


L' = i pUqR jT + i pR J' dx A Fd(^£) - 23 t pi J' RAidx - 

x=l ° ° 

ppR y^s(Ric - - p^^oi ^ (s£c)^ (52) 
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In the case of configurations possessing mirror symmetry about both 
the y and z axes (no dihedral), the general expression of equation (5I) 
can be greatly simplified. If one writes for such a case 

q> = vq>2 + WCP3 + pCp^ + 

where and 9s represent pure translations along the y and z axes, 

?4 is a pure rotation about the x axis, and 9^ is associated with 
the variation of cross section with x, it can be shown from symmetry 
considerations that 


Hence, 


*/ cpguc = jT ® 


= 0 


R 


jT tp£d£ = pR jT cp^^dC 


and, since ^4 is the potential due to pure rotation about the x axis, 
it follows that 

Therefore, recalling that 0 for configurations having symmetry 

about the y and z axes, one finds that the general expression of eq.ua- 
tion (51) reduces to 


L' 


pUopR jT 94Sd5 - 23tpl j^ RAidx 
x=l ° 


(53) 


and it has already been pointed out that for these cases Ai is inde- 
pendent of the rolling velocity p. Thus, for symmetrical configura- 
tions, the rolling moment has now been expressed as the sum of two inde- 
pendent parts, one due to pure rolling and one due to piire translation. 


STABILiry DERIVATIVES 


The specific maneuvers to be considered here will now be defined 
so that stability derivatives can be determined. Rotations are performed 
about the . xyz (body) axes^ at an attitude defined by the angles of 


^Rolling about the wind axis can be treated as a special case by 
the proper choice of the arbitrary body E^xes. 
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attack and sideslip at t = 0, and 
accelerations are permitted in the 
fixed vertical and horizontal 
planes. In particular, if the 
angles of attack and sideslip at 
t = 0 are denoted by Oq and 
respectively, then the velocities 
V and W at any time t are given 
by (see sketch) 

V = UoPoCiOS pt + 

UoCioSin - r(x - ci) 

W = - UoUoCOS pt + 

Uo 3 osin pt - q.(x - Cl) 


e 



From these expressions one finds that 

R = V - IW = (UoPo + iUoOo)e-lPt - (r - ig) (x - ci) 


and 


^ - ip(UoPo + iUoOo)e-ip-t + e"ipt (Uo 3 q + - 


(r - i 4 ) (x - Cl) - ltp(UoPo + iUoOo)e“iP't 

It will be noticed that in the above expressions, the velocity Uq (along 
the X axis) is considered constant. This means that pure pitching and 
yawing motions (q., r O) are performed at constant angles of attack and 
sideslip, so that for such maneuvers the airplane follows a curved fli^t 
path. Now, setting t = 0 in the above expressions, one finds 

R = TJoP t iUoa - (r - iq) (x - ci) 

and 

II = + iUoa - ip(UoP + iUoa) - (r - Iq.) (x - ci) 

and these relations can be substituted directly into equations (^ 5 )^ (^T)j 
and (51) for the forces and moments. It will be noted that the subscript 
on d and 3 has been dropped. This means that for the rolling case 
(p 0) d and 3 of equation ( 5 ^ 4 -) are not the time ^ rates of change of the 
actual angles of attack and sideslip since d and 3 are measured in the 
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fixed vertical and horizontal planes. The alternative would have been 
to define d and 3, as ^ the actual rates of change of a and p (i.e.^ 
a = oq - J3p and p = po + op at t == O) in which case the coupling of 
p with a and p would be obscured by the definition. This would lead, 
for example, to a nonzero value of when the maneuver consists of 

pure rolling at an initial angle of sideslip. This seems undesirable. 

In seeking stability derivatives for the general problem \mder con- 
sideration, it will be advantageous to employ the transformation of the 
airplane cross section to the circle (see Appendix B) . In this way, it 
will be possible to carry out differentiations of the forces and moments 
explicitly with respect to the airplane motions and thus obtain the 
stability derivatives in terms of the transformation without specifying 
the complex potential. Thus, from Appendix B, 


Ai = Rai - Rro^ - Ip ( aoTo^ + a^ao + + S2f|. + 

V 


and therefore 


| 5 . - ^ - lp(at,ro^ + aiaj, + . 

ot ot ot 


where ai is the coefficient of the l/a term of the mapping function 
and ro is the radius of the transformed circle. It_ls recalled that, 
as shown in Appendix B, Ai^ is simply the value of Ai at a = p = p = 

= r = 0 and is therefore directly associated with the shape of the con- 
figuration and the choice of axes. 

The stability derivatives will be obtained by partial differentiation 
of the forces and moments with respect^ to each of the ten independent 
variables a, p, p, g, r, a, p, p, second derivatives will be 

included; that is, there will be derivatives of the types 



• « • 

where all derivatives are evaluated at a=p=p = q = r = a= P = p = 

4 = r = 0. The reason for this choice (which is not customary) will 
become more evident later, but it can be seen at once that all deriva- 
tives defined in this manner are constant for a given configuration and 
that there will be cross derivatives" of the type which will 

show the mutual influence of the longitudinal and lateral motions. Thus 
the total rolling moment due to sideslip, for instance, will be express- 
ible as 
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* _ 


= c 


ip + 


+ 






VCi 


pp 


For the sake of consistency, all coefficients and derivatives are based 
on the same reference area Sr and reference length lr» 


If equations and sere rewritten with the aid of equations 

to (56), the necessary partial differentiations can be effected with 
no difficulty and all the derivatives of the side force, lift, yawing 
moment, and pitching moment are obtained. The derivation of the rolling- 
moment derivatives is not quite so straightforward as the others since 
the expression of equation (51) for the rolling moment contains integrals 


n — 


of the type which cannot be handled directly by residues. 


as mentioned previously. However, when this Integral is differentiated 
with respect to any motion except p or p, the resulting integral can be 
related to one of the integrals already evaluated by residues. In partic 
ular, if we write once more 


9 = vep^ + W93 + pq>4 + cp’ 

it follows from the boundary condition of equation (30) that 


^ = y Z Si 

hn. ds ds 


or 

— rVP 

= y dy + z dz = ds 
on 

Thus the integral appearing in for example, will be 

But, by virtue of Green’s Theorem (see ref. I3, p. 




and again from equation (30) it can be seen that 



26 


NACA TN 3283 


ds 


so that, finally, one can write 


B, f f ‘’’4'^^ = - Uo I f 


This latter Integral has already been evaluated with equations ( 3 ^) and 
(te) to calculate Cy . Similarly, the integrals appearing in the other 

rolling -moment derivatives can also be evaluated and expressed in terms 
of the mapping function of the cross section. The actual differentia- 
tions to obtain the rolling -moment derivatives are simple enough if one 
notes that the order of differentiation of the potential T is important 
in that the expression of equation (51) has been written for a specific 
instant of time (t = O) . That is, T must be differentiated with respect 
to time first, then integrated to give the rolling moment, and finally 
differentiated with respect to the desired motion. Thus, since 9 is 
linear in the angles a and 3 as well as in the angular velocities 
p, q, and r, one observes that 

So8t 8p8t 

and further that 

8oo8t 


8p8t 8(x8p 




5p’ 


= — : etc. 


The resulting expressions for the stability derivatives are given 
in table I, which is arranged so that all the side -force derivatives 
appear in the first column, all the lift derivatives in the second, and 
so on. It is found that a number of derivatives vanish identically, 
that is, regardless of the shape of the cross section. As a matter of 
fact, all but Qk of the possible 325 first and second derivatives vanish 
identically. For obvious reasons, the stability derivatives that are 
identically zero are not listed, but a definite pattern can be seen in 
table I. which shows, for Instance, that all second derivatives of Cl, 
Cy^ Cjn, and vanish except those involving^ p and that there are no 
second derivatives involving a, 3 , p, q, or r. It should be noted that 
the order of differentiation is immaterial so that = ^Lpa^ 
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(if this were not the case, the total number of possible derivatives 
would be 550 . ) In the next section, a total of k$ relationships are 
found among the stability derivatives, so that, at most, only 35 deri- 
vatives need be calculated directly for any given configuration. It is 
important to note that the mapping function must be expanded in exactly 
the form of equation (i^^) before the proper coefficients can be obtained 
for use in the formulas of table I. In particular, the coefficient of 
the cr term must be unity. 


Relationships Among the Stability Derivatives 


Frcan the general results shown In table I, a number of interesting 
reciprocal relationships which are independent of the configuration are 
observed at once. For Instance, the side force due to angle of attack 
is equal to minus the lift due to angle of sideslip * 

Similar equalities among the various derivatives are found throughout 
and are listed in table II. It should be noted at this point that many 
of the relationships . of table II would be obscured by evaluating the deri- 
vatives at a 5^ 0, as is customarily done. 

Beyond these ‘siinple equalities, there are some interesting relation- 
ships which can be brought out by an integration by parts. For example 
(see table l), 

Cj. — r^(x - Cl) A [2itR(ai + ro^) - S] dx 

^ SrlrJo 

= . ^ |(i - ci)[2rtR(i:i + ro^) . [2«R(ai + 

Sr ^ L '-'o ^ 


= - % + 


But from table I it is seen that 



( 


Z - Cx 


) 




Thus, one finds that 
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This is a particularly interesting relationship in that it enables one 
to calculate from the static lift-, and moment -curve slopes a quantity 
which would require dynamic tests in the wind tunnel for a direct exper- 
imental measurement. Another quantity of interest in this category can 
be obtained by Integrating by parts the expression for given in 

table I, The resulting relation is 


- 



(58) 


which states that the damping in pitch is proportional to the lift-curve 
slope. This result, which is independent of configuration or choice of 
axes, was obtained previously by Bryson (ref .12) whose analysis was 
implicitly restricted to bodies having a straight-line axis (i.e., no 
camber of the body or the wings) . The apparent mass concept, which was 
the basis for the analysis of reference 12, will be discussed in a later 
section. Other relations obtained here in a manner similar to that for . 
equations (5?) and (58) are given in table II. 



It is of some interest to look into 
the damping in pitch of wing-body combi- 
nations on the basis of equation (58) . 
Since the lift -curve slope is determined 
entirely by the trailing-edge configura- 
tion (see table l), it is evident that a 
wide variety of airplanes can be treated 
at once quite simply. If one considers, 
for instance, a configuration whose 
trailing-edge cross section consists of 
a circle with symmetrically placed 
straight lines, as shown in the sketch, 
it is known from the transformation (refs. 
k and 12) that the lift-curve slope is 
proportional to the quantity 



Thus, from equation (58), it follows that the ratio of the damping in 
pitch of the wing -body Combination to that of the horizontal wing alone 
is also given by this quantity. It is important to note that any changes 
in shape ahead of the trailing edge (e.g., camber and thickness of the 
wings, variation in fuselage shape, etc.) are immaterial. Th\is, for such 
configurations, the damping in pitch is plotted in figure 1 and it can 
be seen that (l) the body is always, destabilizing, and that (2) this 
effect is a maximum for a body diameter to wing span ratio of 
The damping in pitch is made a maximum, on the other hand, by bringing 
the body to a point (or a line) at or ahead of the wing trailing edge. 



NACA TN 3283 


29 


The relationships given in table II are all Independent of the con- 
figuration, at least to the order of the present analysis. It can be 
seen from table I that for configurations involving any symmetries (for 
example if S5c = there will be additional relationships among the 
stability derivatives. Such a case -will be considered in the section 
titled "Applications of the Theory." 


APPARENT MASS 


Perhaps the most striking featin-e of the results presented in table 
I is the frequent appearance of the quantities [23tR(ai - + S], 

[2jcH,(ai + ro^)* “ S], and 23rl(ai) . In fact, since both ro^ ajid S are 
real, the first two quantities are simply the real parts of 
[2ir(ai - ro^) + S] and [23t(a|L + Tq^) - S] , respectively, while the third 
can be written as the imaginary part of either of these bracketed quan- 
tities. It now becomes evident that a large number of the stability 
derivatives depend only on these two bracketed quantities which, in turn, 
depend only on the shape and size of the airplane cross section. As a 
matter of fact, it can be shown (ref. 17) that the quantities 
pR[2jc(ax - To^y + S] and pR[2jr(ai + r©^) - S] are identical with the 
integrals defined in incompressible flow theory as the additional appar- 
ent mass of the cross section in the y and z directions, respectively. 


These are 


given by jT p9^dz and jTp^gdy (refs. I 3 and 1?) where 9^ 


and 9 q are the velocity potentials for unit velocity of the cross sec- 
tion in the y and z directions. 


The mathematical basis for the use of the apparent mass to calculate 
the transverse force derivatives of slender bodies in steady flow (as done 
by Munk and Jones in refs. 1 and 2) wsls established with Ward’s formula- 
tion (ref. 5 ) of the general expression 


Y 


IL = ipUo jT <Pd£ 
x=Z 


(59) 


since a differentiation with respect to angle of attack, for example, 
yields 


h(Y - IL) 
8a 



Hence, taking the imaginary part of both sides, one finds 
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Itt = pUo^ J <Pady = Uo®(m3a)x=i (60) 

where 11133 Is the additional apparent mass of the cross section in the 
z direction for pure translation in the z direction. It thus becomes 
clear that the lift-curve slope is^ in all cases, given by the apparent 
mass mgg evaluated at the base of the body. If equation (59) "were 
differentiated with respect to angle of sideslip rather than angle of 
attack, then the derivative L3 would be found to Involve rather 

than mQ3, where mgs is the apparent mass in the z direction for pure 
translation in the y direction. 


There remains now the question of the relation between the apparent 
masses and the total forces given by equation (59) ; that is, under what 
conditions can the total forces be calculated from the apparent masses? 
This can perhaps best be clarified by setting 


cp = Vc^2 + WcTg + T ' 


where the first two terms represent rigid-body translations of the cross 
section in the y and z directions and T’ represents variations in the 
shape and size of the cross section with x. For the steady case, the 
velocity components of the centroid of the cross section are given by 




and 


Wc = - Uo 


a - 


dzc 

dx 


so that equation (59) can be written in the form 
Y - ii 


= lUo^ ^3 + - £pE.^(ms3 “ iiQsa) j 


ipUo / 
x=Z 

Note that the quantities 




and 


a - 


<3.Zc 

dx 


represent the angles 


that the line of centroids of the cross sections makes with the flight 
direction (i.e., the local angles of attack and sideslip). It can now 
be seen that the complex force of equation (59) le given by the apparent 
masses and the angles of attack and sideslip of the base cross section, 

provided that jTcp'd^ = 0. 
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Ab an example, consider first a cambered body with no wings whose 
cross sections near the base are adl circular. In this case 


'P'x.I = I' 


T-riJ{y - yc)‘ 


+ (z - Zq)^ so that 6 ^P’d£ =0 and the 


apparent mass gives the side force and lift as well as their derivatives. 
On the other hand, consider a flat plate having deflected partial-span 
flaps. Here (when the plate is alined with the flight direction), the 
entire potential at the trailing edge is given by whose Integral 



does not vanish. 


Hence, in this case, the lift -curve slope is 


given by the apparent mass but the lift is not. This is also true for a 
body of revolution having flat-plate wings at incidence to the body. 


For the unsteady case, as a result of a recently published report 
by Miles (ref. ll), one can show that the stability derivatives can also 
be obtained from apparent mass considerations. This is most easily seen 
from the general expression of reference 11 


Y - U. = IptTo jT <PdS + ip ^y^.dx jlT £pd£ 
x=Z ^ 


since, for example. 


( 61 ) 


^(Y - IL) 





so that the imaginary part yields 

XSr? P ^ 

L4 = -^ / msadx ( 62 ) 

Jo 

It should be noted here again that the total forces themselves are not, 
in general, given by the apparent masses. Miles also shows in reference 
11 that the rolling moment about the wind axis is given by 


L' =|pUottjir Fd(££) +1 play^dx jT Fd(££) (63) 

x=Z ° 

Now, by reasoning exactly parallel to that for the steady case, it can 
be concluded from equations (l2), (6l), and ( 63 ) that in all cases all 
of the stability derivatives (except drag) for rigid slender bodies can 
be obtained frcxn the apparent masses (or, more generally, the "inertia 
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coefficients”). Also obtainable from the inertia coefficients of the 
cross sections are (l) the total side fdrce and lift, provided that 



0, (2) the total yawing and pitching moments, provided that 



0, and (3) the total rolling moment, provided 


that R, 


jT 

x=l 


0. It is interesting to note that conditions (l) and 


(3) Involve symmetries of the cross sections only near the base, but that 
condition (2) is more stringent and is satisfied, for example, by having 
bilateral symmetry of the cross sections over the entire length of the 
airplane . 


An alteinative form for the integrals representing the apparent 
masses or inertia coefficients can be shown to be identical with the 
Integrals representing the kinetic energy of the fluid associated with 
any desired unit velocity (linear or angular) of the cross section; for 


example, jT 


However, it is essential to note 


that it is only for rigid-body motions of the cross section (as repre- 
sented, for example, by cp^) that the two Integrals are identical, since 


only for such motions does 


^n 


- S2L at the boundary. Thus, it can be 
ds 


seen 


that jT p9‘dy jA - pCp’ ds 
dimensions by dividing by a velocity. 


even if is given proper 


Inasmuch as relatively few inertia coefficients have previously been 
calculated, there seems to be little advantage ( other than brevity) in 
expressing the stability derivatives in terms of these coefficients. It 
is felt that the formulas of table I Involving the mapping function will, 
in general, be found more useful, although one should certainly make use 
of any of the coefficients already calculated. In this connection, the 
reader is referred to a recent paper of Kuertl, McFadden, and Shanks 
(ref. 17J in which the apparent masses of a number of interesting cross 
sections are listed for simple translation in the y and z directions. 


The apparent mass integral 


cPsdy was also calculated for a few 


shapes in connection with minimum drag problems in reference I8. How- 
ever, the integrations there were (for the rectangle) carried over both 
the exterior and the interior of the cross section since the configura- 
tion treated there was Indeed a hollow rectangle made up of fo\or thin 
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wings. _ It is necessary then, in order to obtain the desired quantity 
[ 2 irR(ai + r^) - S], to subtract the cross-sectional area (as pointed 
out in ref. 17). It should be pointed out in this regard that all the 
results of the present paper are for solid bodies. Thus, if the body con- 
tains a Jet, one must add to the calculated lift the negative rate of 
change of momentum (in the z direction) of the air passing through the 
Jet. For a simple bent, circular, thin-walled pipe flying approximately 
along its axis, a Jet exit velocity equal to the flight velocity Just 
doubles the lift given for the solid circular cross section. 


The tabulated values of 



in reference 18 have been adjusted 


for this internal flow and additional values have been calculated here to 
extend the range of the variable. The results are plotted in figure 2 to 
show the increase in lift -curve slope and damping in pitch obtained by 
use of a blunt trailing edge of rectangular cross section, of vertical 
end plates near the trailing edge, and of a biplane with sharp trailing 
edges. It will be recalled that both the lift-curve slope and the damp- 
ing in pitch depend only on the traillng-edge cross section, so the 
results of figure 2 are independent of wing thickness, camber, body shape, 
etc., ahead of the trailing edge. 


n fins 

/7 JbJ 


It might be mentioned 
here that if the apparent 
masses of a given cross sec- 
tion in two orthogonal 
directions are equal, then 
the apparent mass of the 
cross section is independent 
of its direction of trans- 
lation. This follows from 
the fact that the momentum 
vectors and the velocity 
vectors add in exactly the 
same fashion. In reference 
17 it was shown that the 
cross sections in the sketch 

possess this Important property. It also follows, then, that ax = 0 so 
that many of the stability derivatives vanish for such configurations 
(see table l) . 




n comers 
regular polygon 


APPLICATIONS OF THE THEORY 


In this section, the results of the foregoing analysis will be 
applied to the calculation of the stability derivatives for several more 
or less special configurations. The first problems to be treated here 
will be concerned with the Introduction of wing thickness and camber as 
parameters since the present analysis is applicable to unsymmetrical 
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configurations. Then, although the finned body of revolution has been 
treated by many authors (e.g., ref. k, 9, 10, and 12), a plane wing-body 
combination will be considered in order to investigate the effect of the 
squared terms in the pressure relation. Finally, stability derivatives 
will be calculated for a wing -body -vertical -fin combination. 


Wings with Thickness and Camber 


The quantity "camber," as Introduced here, may be complex, the 
imaginary part corresponding to the conventional camber in the vertical 
(xz) plane and the real part corresponding to a lateral camber or 
"wiggles" in the horizontal (xy) plane. Perhaps the simplest configura- 
tion of Interest for the present problem can be made up of elliptic cross 
sections whose eccentricity and position in the lateral planes are arbi- 
trary functions of x. The required transformation for such a configura- 
tion is (see sketch) 



C 



ii-cy 

so that 


. = 0 


These quantities 

suffice for the calculation of many of the stability derivatives directly 
from table I. However, for the rolling derivatives, the complex potential 
is, in general, required. The complex potential in the transformed cir- 
cle plane can be derived from reference 16, page 239^ is given by 


®*o “ 


a2 - b^ 


a^ — a« — 


and 


a+b 


the radius of the transformed circle is "2 
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Fi(cr) =: B Zn O’ 


H^»(a + b)(b cos 7 + xe. sin 7 ) 
2 a 


ip(a^ ~ b^)(a + b)^ ^ 

l6o2 A, 

n=e 


( 65 ) 


where and 7 define the lateral velocity of the centroid given by 


Hceiy = R + Uo 



iP^c 


( 66 ) 


The logarithmic term of equation (65) permits a variation in size of the 
ellipse with x. The series term represents the portion of the potential 
required to meet the boundary condition of equation ( 2 l) when dv/dx is 
arbitrary (warped body); this permits a variation of the eccentricity 
b/a with X. If now the potential of equation (65) is transformed to the 
£ plane by means of equation (64) and the coefficient of the l/£ term 
is evaluated at R-p= 0 , itis found that 




2 



— + b — 
dx dx 


£c 




(67) 


With this result, one can obtain all of the stability derivatives except 
and Ci^ directly from the formulas given in tables I and II, for 

any given configuration in this category. Two examples will now be 
considered. 


The first example will consist of a cambered elliptic cone; that is, 
an elliptic body with constant eccentricity whose span is a linear func- 
tion of X. The body axis will be chosen to pass through the center of 
the base and the camber line will be represented by a sine curve 
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sec. A- A 


plotted In figure 3. It Is seen that 


Zc = 6 sin jrx/c as shown In the 
sketch. This body can be thought 
of as a special case of a cam- 
bered triangular wing having a 
blunt trailing edge. The camber^ 
of course, must be sufficiently 
small to Insure the validity of 
the slenderness approximations. 

On the other hand, the thickness 
ratio as Indicated by b/a is 
arbitrary. Thus the range of 
b/a from the flat plate 
(b/a = 0) to the circular cone 
(b/a = 1) can be treated as one 
problem. It should be mentioned 
that the choice of body axes is 
arbitrary so that maneuvers about 
any other set of orthogonal axes 
fixed in the body could be 
handled equally well. A few of 
the interesting stability deri- 
vatives have been calculated for 
this configuration by the formu- 
las of tables I and II and the 
derivatives Cnp and are 

Cnp is always negative for posi- 


tive camber, is linear in the camber, and increases with the thickness. 
On the other hand, Cnap; which is negative in all cases, is independent 

of the camber and decreases to zero as the thickness ratio increases to 
one. 


As a second example of a configuration having elliptic cross sec- 
tions, the "wing-like" shape developed by Squire (ref. 19) is chosen. 

This shape has a variation of eccentricity of the ellipse such that all 
profile sections (except the midspan section) have a rounded leading edge 
and a pointed trailing edge. In addition to this particular thickness 
distribution, for the present problem the wing will be given a camber 
identical with that taken for the elliptic cone; that is, & sin jrx/c. 
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Thus, the wing to be treated 
here has the shape illustrated 
in the sketch. For this wing, 
the major and minor axes of 
the ellipse are given by 

a = ^ X (68) 

and 



where A is aspect ratio and 
t is the maximum thickness 
(thickness of the midspan sec- 
tion at X = c/2) , Some of 
the stability derivatives have 
been calculated for this wing 
by the formulas of tables I 
and II and Cnp and Cn^p are 

plotted in figure 4 . It is 
found that Cnp displays, in 

general, the same variations 
with camber and thickness as 
did that for the elliptic cone. 



sec. A -A 


However, Cn^p displays a trend opposite 


to that for the elliptic cone. That is, Cn^p is seen to Increase with 

thickness, so that the angle -of -attack contribution to the yawing moment 
due to rolling is apparently heavily influenced by whether the trailing 
edge is blunt or sharp. The derivative Cyp was also calculated for 

both wings and was found to be independent of the thickness. In fact, 
for either wing, for the axes chosen. 



While this result appears to contradict the corresponding relation found 
in reference 12 (Cyp = Cxp) ^ it simply highlights the fact that the analy- 
sis of reference 12 does not Include camber although it could be extended 
to do so. Clearly, if 8=0 the two results are in agreement. The 
derivative C^^ was also calculated for both cambered wings and was 

similarly found to be Independent of the thickness. In fact, C^^ was 
found to have a value ecjual to Cnp for zero thickness; that is, 

C^^ = - 0.4508 8 /span. 

It is interesting to note that since the trailing-edge cross section 
of the Squire wing is a straight line, any stability derivatives that 
depend only on the mapping function of the trailing-edge section (see 
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etc.). It is important also tliat for an elliptic cross section. 


the 


quantity [ 2jt(ai + Tq^) - S] (apparent mass in the z direction) is inde- 
pendent of the eccentricity, depending only on the semimajor sixis a, 
while the quantity [2:r(ai - + S] (apparent mass in the y direction) 

depends only on b. Thus it can be seen from table I that all of the 
purely longitudinal derivatives (Ci^, Cl^^, Cl^, Cl^, Cing^, 0^^, 0^^) 

are the same as for a flat-plate wing of local semispan a. Similarly, 
all of the purely lateral derivatives (Cy^, Cn|^^ etc.) are the same as^ 

for a flat -plate wing of semispan b if we replace 3 by minus a, Y 
by L, and N by M. One can see, therefore, that the ellipse is a very 
special cross section and tends to obscure some of the effects of thick- 
ness. For instance, it was seen in figure 2 that [2:tR(ai + ro^) - S] 
for a rectangle Increases with the height of the rectangle and that a 
blunt trailing edge of rectangular cross section will therefore give an 
increase in lift-curve slope over a flat plate and a corresponding 
Increase in the damping in pitch. 


For the evaluation of rolling moments it would appear from equations 
(51) and (52) that some of the integrations might be quite difficult 
because of the nonanalytic character of the integrands. In fact the 
stability derivatives C^p and C^p contain the same nonanalytic inte- 
grands (see table l) . However, the Integrations can sometimes be advan- 
tageously carried out in the transformed circle plane by the method of 
residues. For a configuration having an elliptical cross section at the 
trailing edge, for example, the calculation of_the damping in roll C^p 

becomes quite simple with this technique. Specifically, from table I, 



1 


x=Z 


1 

Sr^r^ 


Rf ^^dcr 
J ^ dcr 


(69) 


Now, since equations (6^1') and (65) for the required transformation and 
complex potential are already in the form of power series in a, and 
since on the circle boundary aa = one chh immediately write the 

integrand as a power series in 0 and therefore use the method of resi- 
dues. It is found for this case that the residue of this series is 


simply i (a<o^ “ ^o^)^ where the subscript refers to the trailing edge. 

8 

Thus the damping in roll is given by 




^r^r^ 



(70) 
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It ie seen from equation (70) that for the Squire wing (bo = O) the damp- 
ing in roll is independent of the wing thickness and camber, remaining 
always at the flat-plate value as given by Ribner (ref. 3)* For the 
elliptic cone, on the other hand, C^p depends on the thickness and the 

effect of a blunt trailing edge beccmies apparent. Values of for 

these cases are plotted in figure 5- should be noted in general that 
since ^F/8p is independent of any translational velocity in the cross 
plane, the damping in roll is not affected by camber. 

The device used above to permit the use of the method of residues 
can always be employed; that is, since ^F/8p, £, and 5 are all express- 
ible as power series in a and a and since a and F are related by the 
radius of the transformed circle, the integrand becomes an analytic func- 
tion of the variable of integration in the transformed plane. However, 
if the transformation itself is an infinite series (as is the case for 
the finned body of revolution), then the residue and consequently the 
damping in roll will emerge as an Infinite series involving combinations 
of all the coefficients of the transformation. This series is, in gen- 
eral, considerably more complicated than that entering into derivatives 
like CYp. 


Plane Wing-Body Combination 


It will be of some interest to consider here certain aspects of the 
plane wing-body problem in view of the fact that some stability deriva- 
tives had been calculated (see ref. lO) before it was generally realized 
that the squared terms in the pressure relation must be retained. A 
number of the simpler derivatives can be obtained quickly from table I 
if the mapping function of the cross section is known, and one of these 
will now be compared with the corresponding derivative obtained in ref- 
erence 10 without the squared terms in the pressure relation. The 
required mapping function is an infinite series obtained by making two 
successive Joukowsky transfomiations (see ref. 4 ), and it is found that 

a© = £c 

"I (• * T-) ■ *“ 

a^ " — a0 ~ . . . 



= 0 
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where s is the wing 
semi span and a is 
the body radius (see 
sketch) . Hence ^ many 
of the stability deri- 
vatives vanish by vir- 
tue of the fact that 
ai is real (see table 
l). It is noted that 
8 , and are here 
ail arbitrary func- 
tions of X (within 
the slenderness 
approximation) and 

the body will have a straight axis (plane wing -body) only if Sc is a 
linear function of x. Also, if Sc is not a linear function of x 
then the evaluation, of Ai^ (which is required for some derivatives) 

becomes a problem which amounts to determining Ai of the complex poten- 
tial for simple translation. This will be done shortly. 








1 




1 





For the pur- 
pose of illustrat- 
ing the influence 
of the squared 
terms in the pres- 
sure relation, the 
rolling moment due 
to sideslip will be 
calculated for the 
special case of a 
flat triangular 
wing mounted sym- 
metrically on a 
cone -cylinder as 
shown in the 
sketch. For this 


case, with the axes chosen as shown, it is clear that ic = 0 , 
a = const. = ao, and s = So(x/c) . It is further noted that (due to 
symmetry) there is no rolling moment provided by the portion of the body 
ahead of x == Xq and that the rolling moment due to sideslip at zero 
angle of attack is zero (C^p = O) . Therefore, the only pertinent deri- 
vative ^ to be calculated is which is given by (see table I) 


^For this configuration all of the coupled (second) derivatives of 
the rolling moment vanish except and ^Zpq^ assumed here 

(for purposes of comparison with ref. lO) that q = 0 . 
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vhere and Zr have been taken as the gross wing area and the maxinnjm 
wing span, respectively. The ratio of this to that for the wing 

alone (ao/so = O) is plotted in figure 6 and is compared with the corre- 
sponding curve of reference 10 . It can be seen that the error incurred 
by the omission of the squared terms in the pressure relation is in excess 
of 100 percent for ratios of body diameter to wing span greater than 0.5. 
Now since = - <^7,^ according to table II, the side force due to 

rolling can also be compared with that obtained from a linear pressure 
relation by means of equation (72). This comparison is presented in 
figure 7 aJ^d it is seen that the difference is even more pronounced than 
that far the rolling moment due to sideslip. It should be mentioned that 
if the wing-body combination is cambered, the contribution of rolling 
moment by the nose will not, in general, vanish, nor will the rolling 
moment at zero angle of attack. 


It has been 
mentioned that for 
a general wing- 
body combination 
of the cross sec- 
tion discussed 
above (i.e., 
not a linear func- 
tion of x), the 
coefficient of the 
l/£ term in the 
complex potential 
must be determined 
if all of the stability derivatives are to be calculated. The coeffi- 
cient Ai for simple translation (p = O) can perhaps most readily be 
obtained from the complex potential for a flat plate in a uniform stream 
by use of the transformation 




b " be 

(see sketch). The complex potential in the w plane is given by (see 
e.g., ref. 16, p. I61) 


Fi(w) = - Ee(w cos 7 - i sin 7 - d^) 


( 73 ) 
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BO that,' by using the above transformation and adding a tenn to 

remove the free stream in the ^ plane and a term B 2n( C ~ £c) allow 
a variation of the radius a with x, one finds after expansion in series 
that the coefficient of l/5 is (for p = O) 

Ai = - Eq j^a^cos y ^ - d^)sln 7 J - b£c 

where d = s + a^/s and it is recalled that EqC^^ = R + Uo(d£Q/dx) for 
pure translation. Now, by noting that for R = 0 


and 


cos 




He V dx/ 

one can write finally for this type of configuration 
Aio = - 


% - r (2a^ - 4^) ^ + a (ye + iZe) 


dx 


dx 


dx 


(75) 


With Aiq determined, all of the stability derivatives except and 

can be obtained directly from tables I and II. The difficulties in 
determining and have been discussed in the preceding section. 


Wing -Body-Vertl cal -Fin Combination 



iT 

be exactly of the form 
reference. 12 is not of 




In reference 12, the mapping 
function was developed for a body 
of revolution having foin* flati- 
plate fins mounted 90^ apart (see 
sketch) . Therefore, one can use 
the formulas of table I directly 
for such a configuration by first 
determining the proper coefficients 
in the expansion of the mapping 
function. It is Important to note 
0 - that for this purpose, the expan- 

sion of the mapping function must 
of equation (U6). The mapping function given in 
this form (as can be verified by carrying out the 
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expansion)^ but it can be modified to give the proper form. The resulting 
transformation is 


£-£)■- . ri (h . I) 


h + f 


where 








and it is found after a somewhat laborious expansion that the first six 
coefficients are 

B-o =~(ja. - f) 


ai = a'^ - 


h + f 


a«= - i(h -f)(^a^ + — 


a« = - a^ + 


,2 . d^ 


h - f 
2 


h + f 

h 


. i (d- + 2 a=)" - ^ 


a (d= + 2a2)^ - I a^j- 






i (d= + 2 a 2 )® + 1 d^a^ + 1 a® + i f 5 -^ Y 

16 U 2 Z \ 2 J 
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In general, 
the Integrations 
indicated in 
table I would 
have to be per- 
formed numeri- 
cally for a given 
configuration of 
the type consid- 
ered above. How- 
ever, for a con- 
ical configura- 
tion, great 
simplification is 
achieved and the 
integrations 
become trivial. 
Ther e fore , a two - 
parameter conical 

configuration is chosen to Illustrate the effect of a vertical fin on a 
number of the stability derivatives; that is, the lower fin is removed by 
setting ts = a, and ti (or t) , a, and s are taken to be proportional 
to X (see sketch). Thus, the coefficients given above and the radius 
of the transformed circle become proportional to x or and most of 
the stability derivatives can be readily calculated. It is noted that 
for this configuration, with the axes chosen,” A^q = = 0 . Also, due to 

the conical property, Cnp = Cm^ = 0 if ci is chosen as 2/3 c. 

The stability derivatives obtainable from the coefficient ai (which 
in this case is real) and from the radius of the transformed circle r© 
have been calculated for a range of the parameters a/s and t/s and are 
plotted in figure 8, The purely longitudinal derivatives are of course 
unaffected by the vertical fin as seen in figure 8(b) (which is inci- 
dentally the same curve as given in figure 1 for a different purpose with 
regard to a more general configuration) . 




Some interesting and important 
effects of the vertical fin can be 
seen in figures 8(a) and 8(c) which 
show a number of the lateral and 
"coupled” derivatives. First it 
should be noted that the two ends 
of the a/s scale correspond to 
the extreme configurations shown 
in the sketch. It is clear that 


a/s^LO the solid curves of figures 8(a) 
and 8(c) have no meeinlng for 

a/s >t/s since this would correspond to a vertical fin inside the body. 
Therefore, for values of t/s less than 1 , the envelope curves (the 
dashed curves) corresponding to a/s = t/s have been plotted to fill in 
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the rest of the range of a/s from 0 to 1. It is seen that even at 
a/s = 0 a vertical fin of small span is quite ineffective on both stabil- 
ity derivatives I and HI (see figs. 8(a) and 8(c)). Hence, the ineffec- 
tiveness of small fins is not entirely caused by the body absorbing the 
fin, but is partly due to a blanketing effect of the horizontal wing. 

The significance of the effect of the vertical fin can perhaps best 
be appreciated by considering a' relatively familiar derivative Cn^^p 

(fig. 8(a)) which can be thou^t of as the angle -of -attack contribution 
to the yawing moment due to rolling. It appears that for fin heights 
such that t/s > 1.0 this derivative becomes positive over some portion 
of the a/s range (almost the entire range if t/s >1.5)* Thus, for 
this configuration, the angle -of -attack contribution to the yawing moment 
due to rolling can be changed from adverse (negative) to favorable 
(positive) by increasing the height of the vertical fin above about I .5 
times the wing semispan, depending on the body radius. 

From the coefficients given in equation (TT)^ it is difficult to 
make any statement regarding the convergence of the series required for 
derivatives such as Cnp (see table l) . Therefore, one is at a loss to 

say how many terms of the expansion must be retained for satisfactory 
calculations. In order to get some idea of the convergence, the deri- 
vatives involving the series were calculated using 2, 3^ 5^ an<i 6 

terms in the series and the results are presented in figure 9* appears 
that for some cases four terms would be sufficient. The results of fig- 
ure 9 indicate a strong blanketing effect of the wing on the vertical fin 
for t/s < 1. It can be seen that for a/t = 1 (plane wing-body) 

Cnp = Cnp = • • . = 0. 

In reference 20, one of the configurations treated corresponds to 
the present conical wing -body-vertical-fin combination for a = 0 (no 
body) . As this furnishes an interesting check on the present calcula- 
tions, the appropriate values of Cyp^ C^p have been taken from 

that report® and are plotted on figures 8 and 9- It caii be seen that the 
agreement is excellent, even for those derivatives calculated with only 
a few terms of the infinite series (fig. 9)- It should be noted in fig- 
ure 9 that wherever the best approximation curve (representing 6 terms 
in the series) cannot be seen, it is because the results were essentially 
identical with the previous approximation. 

It is interesting to note that for a cruciform wing-body combination 
(ti = ts = s), according to equations ( 78 ) and ( 7 ?) and succeeding terms, 
one finds that ao = = a^ = &4 “ ae = • • • = 0 so that many of the 

stability derivatives of table I vanish due to the symmetry of such a 
configuration if the axis of symmetry is chosen as the x axis; for 
example, here again Cnp = Cnp = • • . = 0. 

®The values of Cy , Cyo, and C^g were taken from figures 11, 23 ^ 
and 24 of reference 20 since^there ap^ar to be some typographical errors 
in equation (^) of that report. 
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CONCLUDING REMARKS 


A general analysis has been presented for determining the forces 
(except drag) and moments and the stability derivatives of a slowly 
maneuvering slender wing-body combination of arbitrary cross section. 

The results of the general analysis appear as (l) fornmlas for the forces 
and moments in terms of the airplane shape and motions and *(2) formulas 
for the nonzero stability derivatives in terms of the mapping functions 
of the cross sections. 

Stability derivatives of the first and second order have been con- 
sidered so that the interdependence of the longitudinal and lateral 
motions is Included. A number of relationships among the various stabil- 
ity derivatives were found which are independent of the shape of the air- 
plane, so that, at most, only 35 of a total of 325 first and second 
derivatives need be calculated directly. 

In order to bring out these relationships, the stability derivatives 
have been defined somewhat differently from the usual derivatives. For 
example, the usual derivative of rolling moment due to sideslip would, 
in the present analysis, be given by 

Clp* = Czp + + pcjpp + 

All derivatives as defined in this paper, then, are constant for a given 
airplane. Time rates of change of the angles and angular velocities have 
also been Included, althou^ these effects vanish for the particular 
derivative above. 

The use of the apparent mass concept for problems in slender -body 
theory has been discussed in the light of the present analysis and on the 
basis of previous treatments of slender-body problems by momentum methods. 
It is demonstrated that all of the stability derivatives can be calculated 
from the apparent masses (or inertia coefficients), but that the general 
expressions for the total forces and moments involve additional teimis. 

From the results of the general analysis, some of the stability 
derivatives have been calculated for (l) two triangular wings having 
thickness and camber, ( 2 ) a plane wing-body combination, and (3) a wing- 
body-vertlcal-fin combination. These three cases have been used to show, . 
respectively, (l) the effects of camber, thickness, and bliont trailing 
edge, (2) the influence of the squared terms in the pressure relation, 
and, (3) the effect of a vertical fin on the various stability 
derivatives. 

It was found that the effect of thickness on the angle -of -at tack 
contribution to the yawing moment due to rolling .was essentially opposite 
for blunt and sharp trailing edges, but the effect at zero angle of attack 
was similar in both cases. In both cases, the angle -of -at tack 



NACA ™ 3283 


‘^7 


contribution was independent of the camber while the zero -angle-of -attack 
contribution was linear in the camber. 

The damping in pitch was shown to be proportional to the lift-curve 
slope and therefore to depend only on the base cross section. On this 
basis, the lift-curve slopes and damping in pitch were calculated for 
(1) a wing-body-vertlcal-fin ccmbinatlon, (2) a blimt trailing edge of 
rectangular cross section, (3) a sharp trailing edge with end plates, 
and (^ 4 -) a biplane with sharp trailing edges. Sizable increases over the 
flat-plate values are shown in the last three cases. 

The derivatives usually called Gyp and were calculated to be 

more than 100 percent in error if the squared terms in the pressure 
relation are neglected in the case of a plane wing-body combination of 
body diameter to wing span ratio of O.5 or greater. 

A number of stability derivatives were calculated for a conical wing- 
body-vertical -tail combination and the variations with body diameter to 
wing span ratio were plotted for various vertical tail spans. The influ- 
ence of the vertical fin was found to be markedly altered by the blanket- 
ing effect of the wing for small vertical tail spans. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
June 10 , 1954 
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APPENDIX A 

DIFFERENTIATION OF A CONTOUR INTEC21AL 


z 



Consider a differentiation of the 


contour integral 



•with respect 


to X, where the contour is a smooth 
closed path that depends on x, as shown 
in the sketch. The four points on the 
curve, designated a, b, c, and d, 
represent -the maximum positive and 
negative values of y and z, so that 
the function T on the contour between 
a and c can be designated *?^iow-er 
simply 9^), that between c and a, 

Tupper ^u)^ those between b 
and d and between d and b, 9+ and 
respectively. Thus one can •write 



^a(x) 


9jdy 


+ 


"dx 



Vy 


i -2- / cp^dz - i 

-'b(x) 


pb(x) 

/ 9_dz 

'-'d(x) 


so that there are now four line integrals to be differentiated. This can 
be done directly by means of the formula 


ns(x) 

Sx Jr(x) 


f(x,y)dy 


/ — dy + -2 f(s,x) 

'Jr(x) ^ 


f(r,x) 

OX 


which gives 
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si ■ i.M s ♦ S ■ S * 

rW >Mx,yA(.,ril ^ ^ |. , , . 2. ( . 

'^c(x) OX 8x 8x 

r r^(^) ^[x ^ =^ ^ Y,(x,z)j ^ |d .|b ( 

'^'^b(x) ^x ^x J 

J ri.(») W.[„y,(x,z)l ^ 

LJa.(x) °x J 


where z = Zi(x,y) and y = Yi(x,z) are alternative expressions for the 
contour C and must be single-valued between the prescribed limits of 
integration. 

Now, if the velocity potential 9 is single-valued on the contour, 

then 


c(x) 

and 

9 ^. = 9 « at b(x) and d(x) 

so that all four of the additional terms above cancel and one finds, after 
combining like integrals, that 


— f <pa.r = r ^^[x>y^2i(x,y)I ^ ^ r 89 [x,z,Yi(x,z)] 

Jj, 8x 8x 

But the partial derivatives can be rewritten in the form 

^^[x,y,Zi(x,y)3 ^ ^ ^Zi ~| 

8x “Lsx'^dz Sxl^2i(x,y) 

and 

^9[x,z,Yi(x,z)] r ^9 ^9 ^Yi 1 

8x “LSx'^Sy 8xJy=Yi(x,z) 
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Therefore, since 




-‘z=Zi 


hx 


-»y=Yi 


the final result is 





tiy 




by hx 


dz 






OZ OX 



8y 8x 
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APPENDIX B 

THE RESIDUE Ai OF THE CCMPLEX POTENTIAL 


From eguation (31) it followB that 


2 itg = R£ - R? + ip?£ + 2IU0 


riii 

J dx 


ds + 2 iG(x,t) 


Now, if one maps the body contour in the ? plane onto a circle of 
radius Tq, center the origin, in the 0 plane by the transformation 

e = f(a) = a + £ ^ 
n=o 

(where an is in general complex) then the "boundary function" of 
Milne-Thomson (ref. 16, p. 237 ) is obtained. Noting that on the circle 
boundary 00 = r©^ so that 


r = f (a) = -ffjj. an 


n=o 


one can write the boundary function as 


00 ^ 

Q-n 


n=o n=o 


ip 


('a + £ H- y ^ an ') + T(a) + 2iG(x,t) 


n=o 


n=o 


where 1(0) is simply 2 iUo 




ds expressed as a function of 0, 


It has been demonstrated in reference 16 that this boundary function 
can be satisfied by setting the complex potential equal to the part of 
Q containing only the negative powers of 0. Thus one can set 


F = R 




n=i 


n=o 


m >n 
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where %(a) represents the portion of T containing negative powers of 
o. The residue Ai of F can then be expressed as 

*“ p — ^2®-! a«a2 \ 

Ai = Rai - Rro^ + ip ( aoro*^ + aiao + + — ^ +...)+ Ai^ 

\ ^o ^o / 

where Aiq Is the residue of T|^(cj) and it is noted that Ai^ can depend 

only on the shape of the cross section and not on any of the airplane 
motions; that is, Ai^ _must be a function of x alone and is simply the 
value of Ai at R = R = p = 0. Therefore Aiq Is zero for any con- 
figuration possessing an axis of symmetry if that axis is chosen as the 
x axis. 

It can be seen that Ai is Independent of the rolling velocity p 
provided that ao vanishes and that either all the odd n or all the 
even n are absent in the expansion of f(a). This leads to the con- 
siderations of symmetry given in the text. It should be noted that for 
symmetrical shapes ao is the centroid of the cross section. 
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1 ABIE I.- STABILm DERIVATIVIS - Concluded 
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Figure I 



wing - body - verticat ^fin combi i 
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(a) Blunt trailing edge of rectanguior section^ 



Figure 2,- Damping - in - pitch and lift curve slopes for 
several trailing edge cross sections . 
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Figure 2- 
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(a) Yawing moment due to rolling at zero angle 

of attack . 


Figure 3 . - Effect of camber and thickness for blunt - 
trailing- edge triangular wing (elliptic cone) , 
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fdj Angle -of- attack contribution to the yawing 
moment due to rotting or the rolling 
moment due to yawing . 


Figure 3 - Concluded . 
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t/span 


(a) Yawing moment due to rolling at zero angle 
of attack . 


Figure 4.— Effect of camber and thickness for sharps trailing - edge 
triangular wing (Squire wing). 
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Figure 5 . - Effect of thickness on damping in roll 
for blunt - trailing - edge wing having elliptic 
cross section at the trailing edge . 
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Figure 6, - Influence of the squared terms in the 
pressure relation on calculations of the 
rolling moment due to sideslip . 
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Figure 7. - Influence of the squared terms in the 
pressure relation on calculations of the 
side force due to rolling . 



stability Derivative 



Figure 8 . - Stability derivatives for conical wing- body- 
vertical- fin configuration . 
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Figure 6“ Continued . 



stability Derivative 



a/s 

(c) Stability derivative HI, 



Figure Concluded 





